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Abstract 

I study stellar structures i.e. the mass, the radius, the moment of inertia and the 
oblateness parameter at different spin frequencies for strange stars and neutron stars in 
a comparative manner. I also calculate the values of the radii of the marginally stable 
orbits and Keplerian orbital frequencies. By equating kHz QPO frequencies to Keplerian 
orbital frequencies, I find corresponding orbital radii. Knowledge about these parameters 
might be useful in further modeling of the observed features from LMXBs with advanced 
and improved future techniques for observations and data analysis. 
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1 Introduction 



Presently there are a number of efforts to constrain the dense matter Equations of State (EsoS) 
through astronomical observations of compact stars. The usual approach is to determine the 
mass and the radius of the stars with the help of various observational features like gravitational 
redshifts (z) from spectral lines, cooling characteristics, kHz quasi-periodic oscillations (QPO) 
etc (Lattimer & Prakash 2007, Li et al. 1999, Ozel 2006, 2008, Zhang et al. 2007). But these 
methods are not foolproof, e.g. the value of z used in Ozel's (2006) analysis of EXO 0748—676 
can not be reproduced as mentioned by Klahn et al. (2006). Moreover, to constrain EsoS from 
QPO observations, one needs to believe in a particular model of QPO which is again a subject 
of debate. Another alternative method might be the measurement of the moment of inertia 
from the faster component of the double pulsar system PSR J0737-3039 (Lattimer & Schutz 
2005, Bagchi et al. 2008). Some high mass stars hke PSR J1903+0327, EXO 0748-676 etc. 
prefer stiff EoS and some other stars like 4U 1728-34 (Li et al. 1999), EXO 1745-248 (Ozel 
2008), prefer soft EoS. This fact hints to the possibility of existence of both neutron stars 
and strange stars. But even then, I need some constrains as there are a number of EsoS for 
neutron stars and also for strange stars. Until then, it is interesting to compare the stellar 
properties for different EsoS. For sufficiently fast spinning stars, stellar structures depend upon 
the spin frequency (I'spin)- So the study of stellar structures for rotating stars will help in better 
understanding of the characteristics of fast spinning compact stars like LMXBs and millisecond 
pulsars. That is why here 1 study stellar structures with rotations in section 2. In section 3, 
I study different rotational parameters for strange stars and neutron stars. Although stellar 
structures for rotating neutron stars or strange stars have been already studied by a number of 
groups (some of which 1 discuss in section 3), systematic studies of all relevant stellar parameters 
as well as disk parameters were lacking. That is why here 1 report the variation of a number of 
different stellar parameters as well as disk parameters for different values of star's spin frequency 
and mass in section 2 and section 3. Also I use one EoS for strange stars and another EoS for 
neutron stars whereas in the earlier works people discussed either only neutron star rotations 
or only strange star rotations, there was no comparison between the properties of rotating 
neutron stars and rotating strange stars. In addition I compare my results obtained by using 
a pseudo-Newtonian potential with full general relativistic calculations by other people like 
Haensel & Zdunik (1989) and Lattimer & Prakash (2004) and the close matching found implies 
the correctness of my approach and the validity of the pseudo-Newtonian potential. In section 
4 I discuss a possible application of the knowledge of the rotational parameters in modeling 
kHz QPOs. I end with a discussion in section 5. 



2 Stellar structures with rotation 

I use two sample EsoS of the dense matter among the numerous EsoS available in literature, 
one for the strange quark matter (EoS A or SSA, Bagchi et al., 2006) and the other for the 
nuclear matter (EoS APR Akmal, Pandharipande & Ravenhall 1998). To find stellar structures 
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Figure 1: Variation of the mass with the central density for strange stars (upper panel) and 
neutron stars (lower panel). The parameter is the value of Q, in units of 10^ sec~^. The EsoS 
used are EoS A for strange stars and EoS APR for neutron stars. 
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Figure 2: Variation of the mass with the radius for strange stars (upper panel) and neutron 
stars (lower panel). The parameter is the value of Q in units of 10^ sec""*^. The EsoS used are 
EoS A for strange stars and EoS APR for neutron stars. 
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Figure 3: Variation of the moment of inertia with the mass for strange stars (upper panel) and 
neutron stars (lower panel). The parameter is the value of Q in units of 10^ sec"''^. The EsoS 
used are EoS A for strange stars and EoS APR for neutron stars. 
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Figure 4: Variation of the oblatcncss parameter i.e. the ratio of the polar radius to the equatorial 
radius with the central density for strange stars (upper panel) and neutron stars (lower panel). 
The parameter is the value of VL in units of 10^ sec~^. The EsoS used are EoS A for strange 
stars and EoS APR for neutron stars. 
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Figure 5: Variation of a/Rg with the mass for strange stars (upper panel) and neutron stars 
(lower panel). The parameter is the value of Q in units of 10^ sec""*^. The EsoS used are EoS A 
for strange stars and EoS APR for neutron stars. 
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with rotations, I use the RNS code^. Following Komatsu, Eriguchi & Hachisu (1989), this code 
constructs the compact star models by solving stationary, axisymmetric, uniformly rotating 
perfect fluid solutions of the Einstein field equations with tabulated EsoS (supplied by the 
users) . 

The fastest rotating compact star known till date is probably XTE J 1739-285 (Kaaret et al. 
2007) having z/spin = 1122 Hz, although the measurement has not been confirmed later. The 
second fastest one is J1748-2446ad (Hessels et al. 2006) with Ugpin = 716 Hz. In this work, 
I choose the angular frequency (Q = 27ri'spin) as 2000, 3000, 4000 and 5000 sec'^ (which 
correspond to Ugpin as ~ 318 Hz, 477 Hz, 637 Hz and 796 Hz respectively) . All of the fast 
rotating compact stars except XTE Jl 739-285 have Uspin in that range. I have also computed 
non-rotating, spherically symmetric stellar structures by solving TOV equations which are 
sufficient for slow objects like EXO 0748—676 {I'spin = 45 Hz). Throughout this work, I take 
the stellar mass (M) to be always greater than 1.1 Mq as observations usually hint the stellar 
mass to be greater than that value. 

In Fig 1 I plot the mass against the central density (e^) both for strange stars and neutron 
stars. For a fixed value of ec, M increases little bit with the increase of Q. For all of the values 
of r2, M first increases with the increase of €c (^ > 0) and then after a certain value of M 
(Mmax) starts to decrease (|^ < 0). The stars are unstable when ^ < 0. This instabihty 
appears around ec = 4.1 x 10^^ gm cm~^ for strange stars and around ec = 2.8 x 10^^ gm cm~^ 
for neutron stars; these values does not change more than 5% with the change of Q in the 
chosen range. For any in the chosen range, I get M = 1.1 Mq at ec ~ 1.7 x 10^^ gm cm^^ 
for strange stars and at ~ 0.80 x 10"*^^ gm cm~^ for neutron stars. 

Fig 2 shows the mass-radius plots. With the increase of Q, for both strange stars and 
neutron stars, the maximum mass {Mmax) increases and for any fixed mass, the radius also 
increases due to the larger value of the centrifugal force. Note that here "radius" means the 
equatorial radius R^g which is always greater than the polar radius Rp. For a fixed Q, the 
maximum mass for a neutron star is greater than that of a strange star. For fixed values of Q 
and M, R^q is larger for a neutron star than that of a strange star. The compactness factor 
(M/R) of strange stars is larger than that of neutron stars and the variation of M with R^q 
follows an approximate Rl^ law for strange stars in contrast to neutron star's approximate R~g 
variation. 

Fig 3 shows the variation of the moment of inertia (/) with the mass. For any fixed mass, 
the moment of inertia increases with the increase of Q both for strange stars and neutron stars. 
For fixed values of M and Q, a neutron star possess much higher value of / than a strange star 
because of its larger value of Rgq. 

In Fig 4, I plot the oblateness parameter i.e. the ratio of the polar radius to the equatorial 
radius (Rp/Req) with the mass. It is clear that for a fixed mass, this ratio decreases with the 
increase of fl both for strange stars and neutron stars i.e. the star becomes more and more 
oblate due to the larger values of the centrifugal force. Moreover, for a fixed value of Jl, this 

^http://www.gravity.phys.uwni.edu/rns/ 
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ratio decreases with the decrease of the mass as there the centrifugal force becomes increasingly 
more dominant over the gravitational force. Note that the variation of Rp/ Rgq with M is steeper 
for neutron stars than that for strange stars, but for both of them, the steepness increases with 
the increase of fl. 

In Fig 5, I plot a/Rg with M where Rg = GM/c^ and a = IQ/Mc. As expected from 
their expressions, the plot shows that for any fixed mass, a/Rg increases with the increase of Q 
as expected and for a fixed Q, a/ Rg decreases with the increase of the mass. For the same value 
of M and Q, a neutron star has larger value oi a/Rg than that of a strange star because of its 
larger value of /. a/it!^ is an important parameter of the compact stars as it can be identified 
as the specific angular momentum of the star and its value determines many other properties 
of the star. 

For any other EoS, the value of M^ax and corresponding radius will change depending upon 
the stiffness of that EoS. But the general trend of the M-R curve will remain the same i.e. 
M oc for strange stars and M oc R~^ for neutron stars. The nature of M — ec curve will 
also remain the same. 



3 Rotational parameters 



With the output of the RNS code, i.e. M, R^q and /, I calculate some rotational parameters 
for strange stars and neutron stars. First I calculate the radius of the marginally stable orbit 
which is defined as (Bardeen et al., 1972) : 



where 
and 



Tms = i?5{3+Z2T[(3-Zi)(3 + Zi + 2Z2)]i/2| 

Z, = l + [l-(a/i?,)T/'[(l + a/i?,)^/^ + (l-a/i?,: 
Z, = [?,{a/Rgf-rZlY'^ 



lV3l 



(2) 



(3) 



The "-" sign in the expression of r^s implies the co-rotating motion and the "+" sign implies 
the counter-rotating motion which I call as r^s, co and Tms, counter respectively. As the values 
oi a/ Rg are always very small, both Zl and Z2 have their values ~ 3. 

The Keplerian frequency of a particle orbiting around the star at a radial distance r can be 
expressed as 



1/2 



(4) 



where Fm{r) is the force per unit mass. As an example, I take Fm{r) as derived from a pseudo- 
Newtonian potential by Mukhopadhyay & Misra (2003) 
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Figure 6: Variation of r^s, co with the mass for strange stars (upper panel) and neutron stars 
(lower panel). The values of Q in units of 10^ sec~^ are shown. The EsoS used are EoS A for 
strange stars and EoS APR for neutron stars. 



RgC 



,2 



-(v) + (v)^ 
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In Figs 6 and 7, I plot Vms, co and Vms, counter respectively with the mass. For any EoS, 
T^ms, CO is always smaller than r^s, counter for any fixed values of Q, and M. For a fixed Q, both 
T^ms, CO and Vms, counter increases linearly with the increase of M. But for a fixed M, r„is^ co 
decreases with the increase of Q and r^s, counter increases with the increase of Q. For strange 
stars, both Vms, co and rms, counter are always greater than Rgg. For neutron stars, Vms, co and 
r^ns, counter are smaller than Rgg for low masses and they become greater than Rgg for higher 
masses; rms, co becomes equal to Reg around 1.5 — 1.7 Mq and rms, counter becomes equal to Rgg 
around 1.15 - 1.2Mq. 

All these happen because of the nature of the terms in the expression of r^s (Eqn.l). As 
both Rg{S + Z2} and Rg{[{3 — Zi){3 + Zi + 2Z2)]^^'^} are positive quantities and the second 
one is smaller than the first one, their sum {rms, counter) must be greater than their difference 
{rms, co)- As for a fixed value of Q, both Rg{i + Z2} and Rg{[{3 - Zi){3 + Zi + 2^2)]^/^} 
increases with the increase of M with the first term having much steeper slope, both rms, co 
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Figure 7: Variation of r^s, counter with the mass for strange stars (upper panel) and neutron 
stars (lower panel). The parameter is the value of fl in units of 10^ sec"-*^. The EsoS used are 
EoS A for strange stars and EoS APR for neutron stars. 



and Tms, counter incrcase with increase of M. For a fixed mass, Rg{3 + Z2} remains almost 
constant but -Rc,{[(3 — ^i)(3 + Zi + 2^2)]^^^} increases with the increase of il, so their sum 
(T^ms, counter) increases with the increase of and the difference (r^s, co) decreases with the 
increase of fl. 

For strange stars (EoS A), Rg{3 + Z2} ^ Req for any value of M. So after addition 
or subtraction of a comparatively small term (Rg{[(3 — Zi){S + Zi + 2Z2)Y^^}) with it, the 
expression {vms, counter or Tms, co) remain always greater than Req 

For neutron stars (EoS APR), Rg{3 + Z2} < Req at lower values of M where the values 
for Req arc sufficiently larger, but -Rg{3 + Z2} > Req at larger vahics of M. Here Rg{[{3 — 
Zi)(3 + Zi + 2Z2)Y^'^} is very small in comparison to both Rg{3 + Z2} and Req for all values of 
M. So the addition or subtraction of Rg{[{3 - Zi){3 + Zi + 2^2)]^/^} with Rg{3 + Z2} (to get 
Tms for counter-rotating or co-rotating motions respectively) does not change the overall trend, 
only the addition (for counter-rotation) shifts the transition towards lower values of M whereas 
the subtraction does the reverse thing. On the other hand, for strange stars, Rg{3 -\- Z2} is 
always much greater than Req and even after the addition or subtraction of the smaller term 
Rg{[{3 — Zi){3 + Zi + 2^2)]^/^}, it (rms) remains greater than Rgq. 
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Figure 8: Variation of fjk(r) with the mass for strange stars (upper panel) and neutron stars 
(lower panel) for = 0.4 x 10^ sec~^. Stellar masses are taken as 1.1 Mq and 1.5 Mq. The 
EsoS used are EoS A for strange stars and EoS APR for neutron stars. 
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In Fig 8, I plot the variation of Uk{r) as a function of r. No significant difference in Ukir) 
(at any chosen r) between a strange star and a neutron star having the same values of M and 
fl is observed (specially at higher values of r). In comparison to the whole range of ^'^(r), 
the differences between ^'^(r) for the co- rotating and the counter-rotating motions (keeping 
all of the other parameters fixed) and the variation of i^kif^) with Q (keeping all of the other 
parameters fixed) arc very small. | [uk, co{r) - Vk, counteAr)] / [vk, counter{r)] \ varies around 0.2 
to 0.001 for r = to 500 km depending slightly upon the choice of the EoS, M and 
I [i^fe, oAif) — Vk, o.2(?^)] / [i^k, o.2{r)] I varies 0.1 to 0.001 for r = to 500 km (where the 3rd 
parameter in the subscript denotes the value of fl in units of 10"^ sec~^) depending slightly 
upon the choice of the EoS, M and the direction of the motion. So as an example, the plot in 
Fig 8 is only for co-rotating motion with Q = 0.4 x 10^ sec~^ and for only two chosen values 
of M. 

But these differences between i'k,co{'r) and i^k,counter{f) become larger for smaller values of r 
where i'k,co{r) < i^k, counter (r) , but at sufficiently higher values of r, ffe,co(r) > i^k.counterir). This 
transition occurs at r ~ 26 — 27 km for both strange stars and neutron stars with M = 1.5 Mq 
and Q = 0.4 x 10'^. Similar trends have been noticed for other values of M and Q. The value of 
r at which this transition occurs increases slightly with the increase of M but does not depend 
significantly with Q. 

We will now study the variation of Vk{Req) i-e. t'fe(r) at r = Reg with M taking different 
values of Q and both co-rotating and counter-rotating motions. i'k,co{Req) is the Keplerian 
frequency of a particle orbiting the star at the star's equatorial surface and I obtain it by 
equating the gravitational force with the centrifugal force. 

In Figs 9 and 10 I plot Vk{Req) with the stellar mass for the co-rotating and the counter- 
rotating motions. co{Req) > ^k, counter{Req) always. Let us concentrate mainly on Fig 9 
as Uk, co{Req) can represent the rotational frequency of both the constituent particles of the 
star and the accreting material at the equatorial surface of the star whereas Uk^ counter{Req) can 
represent only the accreting material. I see that Uk, co{Req) > ^spin in the chosen range of M 
and Vspin- Fig- H shows that the difference Au — Uk, co{Req) — ^spin decreases at higher lyspin 
and/or lower M. The condition Au = is the "mass shedding limit" as for Uk^ co{Req) < ^spim 
matter from the star would fiy way due to the centrifugal force. The spin frequency where the 
"mass-shedding" starts (i.e. Au = 0) is the maximum possible rotational frequency of the 
star and it is known as the Keplerian frequency of the star (i/^ = Ugpin = I'k, co{,Req)) and the 
corresponding angular frequency is called as the Keplerian angular frequency {VL],) of the star. 
For fixed values of M and Vspm, AvgSA > Avapr- This imphes that the "mass shedding hmit" 
will be reached in case of neutron stars at comparatively a lower value for Vspm than for strange 
stars which means that strange stars are more stable than neutron stars against rotations (we 
have checked that this fact remains true even if one uses BAG model EoS for strange stars). 
The Keplerian angular frequency (fi^) for compact stars has been studied in the literature by 
different groups like Friedman, Ipser & Parker (1986); Haensel & Zdunik (1989); Lattimer et al. 
1990 and many other authors. 

Haensel & Zdunik (1989) proposed an analytic relation between il^ and maximum allowable 
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Figure 9: Variation of Uk, co{Req) with the mass for strange stars (upper panel) and neutron 
stars (lower panel). The parameter is the value of Vt in units of 10"^ sec"^ corresponding to Vspin 
318 Hz, 477 Hz, 637 Hz and 796 Hz respectively. The EsoS used are EoS A for strange stars 
and EoS APR for neutron stars. 
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Figure 10: Variation of u^^ counter{Req) with the mass for strange stars (upper panel) and neutron 
stars (lower panel). The parameter is the value of f2 in units of 10"^ sec~^ corresponding to fspin 
318 Hz, 477 Hz, 637 Hz and 796 Hz respectively. The EsoS used are EoS A for strange stars 
and EoS APR for neutron stars. 
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Figure 11: Variation of — co{Req) ~ ^spin with the mass for strange stars (upper panel) 
and neutron stars (lower panel). The parameter is the value of VL in units of 10^ sec"''". The 
EsoS used are EoS A for strange stars and EoS APR for neutron stars. 



16 



static mass {Mmax) and the corresponding radius (Rmax) 



( 



M, 



) 



0.5 



( 



) 



-1.5 



ni = 7.7 X 10^ 



Mr 



max 



10 km 



max 



(6) 







With their choice of neutron star EsoS, Lattimer et al. (1990) found Ql to lie between 0.76 x 
10^ — 1.6 X 10^ sec~^ from Eq. (6). Using Bag model EsoS for strange stars, Prakash, Baron & 
Prakash (1990) found ^1 ^ 1.0 x 10^ sec"^ for M > I.UMq. 

Afterwards, Lattimer and Prakash (2004) gave a more useful expression of fl^ of a star of 
mass M and non-rotating radius R. 



To test these two simple analytical expressions (Eqns. 6, 7), I run the task "kepler^^ in RNS 
code which produces stellar configurations for stars rotating with fi^. Using the RNS outputs 
{i.e., M, R, I, Ql, I calculate i/fe, co{Req) from Eqn. 4. At fll, I should get fl^ = STrz/fe, co{Req) 
(by the definition of fll). 

In Fig 12, I plot the variation of with mass. The fine labeled (1) is the maximum limit of 
Ql given by Haensel & Zdunik (1989) i.e. Eqn 6. The curve labeled (2) is from the analytical 
expression given by Lattimer & Prakash (2004) i.e. Eqn 7. The curve labeled (3) is the output 
of the RNS code and the curve labeled (4) is 27ri/fc^ co{Req)- 

For neutron stars (APR EoS) the curves (2), (3), (4) are very close to each other which 
supports the correctness of both the analytical expression of Lattimer & Prakash (2004) (Eqn 7) 
and the Pseudo-Newtonian Potential of Mukhopadhyay & Misra (2003). The maximum value 
of Ql obtained using Eqn. (6) is ~ 11519 sec~^. Here the Kerr parameter a/Rg lies around 
0.66 which is much greater than the values at lower frequencies (Fig. 5) and the neutron stars 
are very oblate having Rp/ R^g ~ 0.59 — 0.56. 

For EoS A, the RNS code failed to perform the task "kepler" for Cc < 1.51 x lO^^gm cm"^ 
and at ec = 1.51 x lO^^gm cm~^, I get M — 1.61 Mq which is greater than M^ax for static 
configuration. Moreover, at this mass, the value of the Kerr parameter a/Rg is greater than 1, 
which is unphysical. At e^- ^ 1.58 x lO^^gm cm~^ (M ^ 1.72 Mq) , a/Rg becomes less than 1, 
but very high (~ 0.9). So no direct comparison of fll obtained with the RNS code with the 
analytical expressions are possible and in the figure I only plot the analytical expressions. Using 
Eqn. (7), I get fll — 13086 sec~^ at €c — 1.51 x lO^^gm cm~^ and increases with the increase of 
Sc (or M). This value of is less than the maximum value of Ql obtained using Eqn. (6) which 
is 14940 scc^^. Here the strange stars are very oblate with RjJ R^q ~ 0.38 — 0.47. To check 
whether all these facts are intrinsic to strange star properties or depend upon the particular 
model of strange stars, I have even checked with the BAG model having model parameters as : 
B — 60.0 MeV/fm^, = 150.0 MeV, = nid — 0, a^ — 0.17 where B is the Bag parameter, 
m^t, m„. 'm,i arc masses of s, u and d quarks respectively giving static Mmax = 1-82 Mq. Here 
also I get a/Rg > 1 at lower ec, a/Rg < 1 when ^ 0.82 x lO^^gm cm'^ (M ^ 2.33M0). Using 
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Figure 12: Variation of Ql with mass. The hne labeled (1) is the maximum limit of Ql given 
by Haensel & Zdunik (1989). The curve labeled (2) is from the analytical expression given 
by Lattimer and Prakash (2004). The curve labeled (3) is the output of RNS code and the 
curve labeled (4) is 27ruK, co{Req)- The EsoS used are EoS A for strange stars and EoS APR 
for neutron stars. 

Eqn. (7), I get fll = 8756 sec^^ at Cc = 0.82 x lO^^gm cm~^ which increases with increase 
of ec (or M). This value of fi^ is less than the maximum value of fll obtained using Eqn. (6) 
which is 10695 sec"^. The strange stars are very oblate having Rp/ Rgq ~ 0.38 — 0.42. 

The probable fastest spin frequency of a neutron star (XTE J1739-285) is Vspin = 1122 Hz 
or Q = 7049.734 sec^^, which is less than the value of Ql of both strange stars and neutron 
stars as derived from Eqn.( 6) or even less than as derived from Eqn. 7 for a canonical value of 
the stellar mass as M = 1.5 Mq (see Fig. 12). So there is no problem of this star being either 
a strange star or a neutron star. One needs to conclude about its nature by other observational 
evidences. 

Several other people performed numerical studies on structures of rapidly rotating compact 
stars like Friedman & Ipscr (1992). Wcbcr & Glcndening (1992); Cook, Shapiro & Teukolsky 
(1994); Eriguchi, Hachisu & Nomoto (1994); Saigado et al. (1994a,b); Gourgoulhon et al. 
(1999); Bhattacharyya, Thampan & Bombaci (2000); Bombaci, Thampan & Datta (2000), 
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Gondek-Rosinska et al. (2000); Bhattacharyya & Ghosh (2005); Hacnsel, Zdunik & Bejger 
(2008); Haensel et al. (2009). In these, people usually kept themselves confined in studying 
the properties of stars hke M, ec, Req, T/W etc or at most properties particles rotating at 
the stellar surface. But the use of the pseudo-Newtonian potential enables us to study the 
properties of particles orbiting around the star at a much higher distance (r >> Rgg) and 
I report quantities like rms,co and r„is,counter, J^k,co{r), i^k, counter (r) which arc useful to study 
accretion onto rotating neutron stars or strange stars (see next section). I also report various 
stellar parameters M, R, I, Rp/R^q, a/Rg, i'k,co{Req), ^k,counter{Req) etc. There is no previous 
work where all this parameters were reported together. Also I use one EoS for strange stars 
and another EoS for neutron stars whereas in the earlier works people discussed either only 
neutron star rotations or only strange star rotations. I also discuss about kHz QPOs within 
the scenario of Mukhopadhyay et al. (2003) in the next section. 

4 Applications in kHz QPO models 

kHz QPOs in the LMXBs are very interesting phenomena. The most popular model to explain 

kHz QPOs is the beat frequency model (Strohmayer et al. 1996) which assumes the upper 
QPO (uup) as the Keplerian orbital frequency z/fc(r) of the innermost orbit in the accretion disk 
around the star, the separation between the two peaks i.e. Az/peafe = ^up — ^low as yspin and 
the lower QPO {uiow) as the beat frequency of u^^r) and Ugpin. This model suggest that for a 
particular star Ai/peafe — ^spin will remain constant. 

However, it has been observed that in many sources, Aup^ak changes with time. This 
phenomenon discards the reliability of the beat frequency model for the QPOs. Two examples 
of alternate models can be found in Titarchuk and Osherovich (1999) and Mukhopadhyay et 
al. (2003). Both of these later two models suggest uigyj = Vkir). But in the present work, I 
don't prefer any particular QPO model over the others. 

We take two sample LMXB for which kHz QPOs have been observed, namely KS 1731—260 
with Pspin = 524 Hz and 4U 1636—53 lygpin = 581 Hz. As none of the above mentioned models for 
kHz QPOs is established beyond doubt, I fit both u^p and uiow to z/fc(r) to find the corresponding 
radial distance r — (from Eq. 4) which I call r^^ and r^^ low respectively. Mukhopadhyay 
et al. (2003) calculated r^, low for these two sources taking mass-radius values for non- rotating 
strange stars and using a/Rg (J in their notation) as an parameter. Here I get rk, up and r^, low 
taking exact values of R^q and a/Rg as generated by the RNS code with the specific values 
of i^spin and chosen values of M (as masses of these two stars are not well determined). QPO 
frequencies for KS 1731—260 has been taken from Wijnands and van der Klis (1997) and those 
for 4U 1636—53 has been taken from Jonker, Mendez & van der Khs (2002). Note that the 
QPO frequencies of 4U 1636—53 shift with time and hence I have taken only one set of values 
as an example. 

Table 1 shows the rotational parameters for KS 1731—260 and 4U 1636—53. Table 2 shows 
the values of the radius of the Keplerian orbit (r^) obtained by equating z/fe(r) with QPO 
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Table 1: Rotational parameters for two LMXBs. 



source 




M 


EoS 




/ 


Rp/ Rgq 


a/Rq 


"^ms, CO 


'^ms, counter 




(Hz) 






(km) 


(lO'^'' gi)i cni^) 






(km) 


(km) 


KS 1731-260 


524 


1.1 


SSA 


7.344 


0.564 


0.982 


0.174 


8.801 


10.656 






1.1 


NS 


11.874 


0.975 


0.918 


0.301 


8.081 


11.297 


4U 1636-53 


581 


1.1 


SSA 


7.348 


0.565 


0.978 


0.194 


8.694 


10.754 






1.1 


NS 


12.015 


0.990 


0.897 


0.339 


7.860 


11.485 


KS 1731-260 


524 


1.5 


SSA 


7.688 


0.908 


0.984 


0.151 


12.179 


14.366 






1.5 


NS 


11.631 


1.191 


0.910 


0.218 


11.135 


15.012 


4U 1636-53 


581 


1 .5 


SSA 


7.699 


0.911 


0.981 


0.168 


12.051 


14.485 






1.5 


NS 


11.718 


1.508 


0.926 


0.278 


11.205 


15.244 


KS 1731-260 


524 






















1.9 


NS 


11.255 


2.032 


0.959 


0.210 


14.894 


18.723 


4U 1636-53 


581 






















1.9 


NS 


11.310 


2.015 


0.919 


0.235 


11.617 


18.9357 



frequencies. I choose the stellar mass to be 1.1 Mq, 1.5 Mq and 1.9 Mq. For a given star 
{i.e. Vspin is fixed), the values of depend significantly on (a) whether I consider the co- 
rotating or the counter-rotating motion, (b) on the choice of the EoS and (c) on the chosen 
value of M. But the values of depend significantly only on M when I keep Vspin and fixed. 

For KS 1731—260, r^^ — fk, up km for both strange stars and neutron stars (indepen- 
dent of the choice of mass) and for both co-rotating and counter-rotating motions whereas for 
4U 1636—53 this value is ~ 4 km. 

f^k, low — Tms ~ 4 to 7 for KS 1731—260 and ~ 7 to 10 km for 4U 1636—53 for co-rotating 
motions, ~ 2 to 9 km for KS 1731—260 and ~ 8 to 10 km for 4U 1636—53 for counter-rotating 
motions; whereas r^^ „p — ~ 1 to 5 km for KS 1731—260 and ~ 2.5 to 6 km for 4U 1636—53 
for co-rotating motions Tk, up — fms ~ —1.5 to 2.5 km for KS 1731—260 and ~ —0.5 to 3.5 km 
for 4U 1636—53 for counter rotating motions. These values decrease with the increase of M 
and note that at M = I.9M0, r^, ^p. counter < f^ms- Remember, Mmax < 1.9M0 for EoS A for 
the values of Vl of KS 1731-260 and 4U 1636-53. 

5 Discussion 

We study the rotational parameters for strange stars and neutron stars using the RNS code 
which consider the generalized axisymmctric metric for rotating stars. I find that the rotational 
parameters like M^^x, Reg, I, Rpl Req, <^/Rg depend on Vspin a-^d the choice of the EoS. The 
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Table 2: Keplerian radius (r^) by fitting kHz QPOs for two LMXBs. 



source 






M 


EoS 


^fc, low, CO 


^fc, low, counter 


'^k, up, CO 


^fc, up, counter 




(Hz) 


(Hz) 


(Me) 




(km) 


(km) 


(km) 


(km) 


KS 1731-260 


898.3 


1158.6 


1.1 


SSA 


15.135 


15.341 


12.923 


13.245 








1.1 


NS 


15.099 


15.442 


12.835 


13.379 


4U 1636-53 


688 


1013 


1.1 


SSA 


18.030 


18.098 


14.018 


14.310 








1.1 


NS 


18.061 


18.166 


13.953 


14.448 


KS 1731-260 


898.3 


1158.6 


1.5 


SSA 


17.071 


17.477 


14.758 


15.255 








1.5 


NS 


16.958 


17.618 


14.609 


15.420 


4U 1636-53 


688 


1013 


1.5 


SSA 


20.080 


20.378 


15.891 


16.396 








1.5 


NS 


20.023 


20.504 


15.748 


16.572 


KS 1731-260 


898.3 


1158.6 




















1.9 


NS 


18.766 


19.631 


16.346 


17.307 


4U 1636-53 


688 


1013 




















1.9 


NS 


21.896 


22.662 


17.507 


18.530 



value of Tms depends on the stellar mass, the choice of the EoS and whether the motion is co- 
rotating or counter-rotating. The dependence of t'fc(r) on the choice of EoS, Q and the direction 
of motion is prominent only at low r Req)- That is why, the values of as obtained by 
fitting the kHz QPO frequencies do not depend much on the choice of the EoS and on the 
direction of motion as here > Rgg. 

The beat frequency model suggests that r^^ is the radius of the innermost Keplerian orbit 
of the {vin) of the accretion disk whereas the models by Titarchuk and Oshcrovich (1999) and 
Mukhopadhyay et al. (2003) suggest that Vk, low = ^in- The disk parameter rj„ can be estimated 
by X-ray spectral analysis. As relativistic broadening is more dominant at the innermost edge 
of the disk, the existence of a relativistically broadened iron Ka line helps one to determine 
the value of rj„ (Reis, Fabian & Young 2009; Di Salvo et al. 2009 and references therein). It 
is clear from the table 2 that for a fixed set of Af and Vspin^ f^in is different for different EsoS, 
but the difference is always < 1%, whereas presently rj„ is measured only accuracy up to 20% 
(Reis, Fabian & Young 2009; Di Salvo et al. 2009). Moreover, one should remember that this 
determination of rj„ depends upon the QPO model which is not beyond doubt at the present 
moment. Indeed there is a strong need of better explanation of kHz QPOs to understand 
various features like the shift of the peaks, their correlations (Bclloni, Mcndcz and Homan 
2007; Yin et al. 2007), side-bands (Jonker, Mendcz and van dcr Klis 2005) etc. Considering 
all these facts, 1 conclude that, to constrain dense matter Eos by measuring Tin using X-ray 
spectral analysis, I need much better accuracy and hopefully future advanced technology like 
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ASTROSAT will provide us such ultra-high accuracy. 

In this work I have primarily confined myself to a maximum value of lygpin as 796 Hz. But 
then I have also studied the properties of the compact stars if they rotate with the maximum 
spin frequency i.e. the Keplerian spin frequency. 
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